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Anti-solar differential rotation and surface flow pattern on UZ Librae 
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We re-investigate UZ Librae spectra obtained at KPNO in 1998 and 2000. From the 1998 data we compose 1 1 consecutive 
Doppler images using the Ca 1-6439, Fe 1-6393 and Fe 1-641 1 lines. Applying the method of average cross-correlation of 
contiguous Doppler images we find anti-solar differential rotation with a surface shear of a —0.03. The pilot application 
of the local correlation tracking technique for the same data qualitatively confirms this result and indicates complex flow 
pattern on the stellar surface. From the cross-correlation of the two available Doppler images in 2000 we also get anti-solar 
differential rotation but with a much weaker shear of a ~ —0.004. 
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1 Introduction 

In our time, the Sun is the only star for which the sur- 
face differential rotation and other surface flows can be 
measured directly. The detection of differential rotation on 
active stars, however, is a challenging observational task. 
Starspots are regarded as manifestations of solar-like active 
regions and as such, are suitable for tracking and recovering 
surface flows. In this paper we make one of the first attempts 
at recovering surface flows on spotted stars (see also the pa- 
per by Kovari et al. 2007a). 

Our target is the RS CVn-type active giant star 
UZ Librae. The star has a very fast rotation rate of Prot — 
4.76 days, which is equivalent to the orbital period. Doppler 
images of the star were presented by Strassmeier (1996) 
and Olah, Strassmeier & Weber (2002). In the latter pa- 
per the signature of a weak anti-solar differential rotation 
(DR) was found but the authors regarded the detection in- 
conclusive. The possible DR was also examined by Olah 
Jurcsik and Strassmeier (2003) using photometric observa- 
tions. They found two permanent, distinct periods in the 
Fourier-spectrum of the long-term data, which were as- 
signed to the different rotation rate of spots at high and 
low latitudes, resulting in a = Afi/fi of —0.0026. Adopted 
stellar parameters of UZ Lib are summarised in Table [T] the 
values which are not marked were taken from Olah, Strass- 
meier & Weber (2002). 



Table 1 Stellar parameters of UZ Librae 



spectral type ^ 


KOIII 


p , - p 2,3 


4.768241 days 


V Sin t 


67 ± 1 kms-i 




50°±10° 


Toff 


4800 K 


log 9 2 


2.5 


Microturbulence ^ 


0.8 km s-l 


Macroturbulence C^r = i^t 


4 km s^^ 


Chemical abundance 


solai' 



* Corresponding author: e-mail: k.vida@elte.hu 



1 Bopp&Stencel(1981) 

2 Strassmeier (1996) 
^Fekel et al. (1999) 



2 Data processing 

We used two sets of spectra obtained at Kitt Peak National 
Observatory (KPNO) in 1998 and 2000. From the 16 spec- 
tra obtained in 1998 we recover a time-series of 11 Doppler 
images. For the first image we take the first six spectra, and 
run our inversion code TempMap (Rice et al. 1989) for the 
Ca 1-6439, Fe 1-6393 and Fe 1-641 1 mapping lines. Then, we 
proceed with the next set of six spectra, leaving the first 
spectrum and adding the subsequent one (i.e, the 7th) to 
the modelled datased, producing a slightly different Doppler 
image. This pattern is followed until we get the last image 
from the last six spectra of the time-series. The respective 
Doppler images for the three different mapping lines are in 
quite good agreement with each other and are averaged for 
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further use, resulting in a time-series of 11 Doppler images 
covering about 20 days (4.2 rotation cycles). Note, that a 
similar method was applied to get a motion picture of the 
surface spots on HR 1099 (Strassmeier & Bartus 2000) and 
later on LQHya (Kovari et al. 2004). From the 17 spectra 
gathered in 2000 we reconstructed only two images from 
two subsets consisting of 8 and 9 spectra, both covering 
about 8 days, but separated by a time lag of 16 days (Olah, 
Strassmeier & Weber 2002). 



3 Results from the cross-correlation study 

The method of 'Average Cross-CORrelation of contigu- 
ous Doppler images' (hereafter ACCORD) uses a combi- 
nation of independent image pairs (i.e., which do not have 
common spectra), whose relevant latitude strips are cross- 
correlated with each other When having multiple cross- 
correlation maps, the elapsed time lags between the cor- 
related images are different, so in order to average them, 
one has to normalize the cross-correlation maps (i.e., ap- 
plying a linear stretching or compression along the longi- 
tude shift) depending on the time lag. Unfortunately, the 
larger the time lag between the cross-correlated images, 
the greater the risks of masking the DR pattern by indi- 
vidual spot evolutions (waxing or waning of spots, emerg- 
ing of new spots, melting, etc.). Thus, to minimize this 
effect, we correlate only the consecutive but contiguous 
(i.e., the closest independent) images. The average cross- 
correlation map is then fitted by a quadratic differential ro- 
tation law: = ricq — Ailsin^ d, where stands for 
latitude, il^q means the angular velocity on the equator, and 
Ar2 = ricq — ^^poic the angular velocity difference between 
the equator and the poles. The differential rotation parame- 
ter is derived as a = Ail/ileq. For a detailed description of 
this method, as well as for further applications see Kovari 
et al. (2004, 2007b). 

The resulting cross-correlation maps are plotted in 
Fig.IUfor 1998 (top) and for 2000 (bottom). The best coiTe- 
lated longitude shifts are fitted by a quadratic DR law. The 
best-fit rotation function indicates zero longitude shift of the 
equatorial zone, i.e., the equator is rotating with the orbital 
period. This result agrees well with the findings in Olah, 
Jurcsik & Strassmeier (2003). The derived a parameters are 
-0.027 ± 0.003 and -0.004 ± 0.001 for the 1998 and 2000 
data, respectively, indicating a weak anti-solar differential 
rotation in both cases, matching (at least qualitatively) the 
negative DR parameter derived from photometric data by 
Olah, Jurcsik & Strassmeier (2003). 

The cross-correlation study can be repeated for the cor- 
responding longitudinal stripes (meridian circles). This case 
we take only the hemisphere of the always visible pole, be- 
cause Doppler imaging is more reliable for this part. We 
then look for the correlation pattern in the same way when 
searching for DR from longitudinal cross-correlation. The 
resulting latitudinal cross-correlation pattern is shown in 
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Fig. 1 Cross-correlation maps obtained from 1 1 time- 
series Doppler images in 1998 (top) and from 2 Doppler im- 
ages in 2000 (bottom). The better the correlation the darker 
the shade. Cross-correlation functions per latitude are fit- 
ted by Gaussian curves (Gaussian peaks are indicated by 
dots, the FWHMs by error bars). The continuous line shows 
the fitted quadratic differential rotation law. Longitude shifts 
from the LCT method are overplotted as triangles in the top 
panel (see Sect.|4]for details). 

Fig. [3] (top). The best correlated latitude shifts indicate in- 
coherent flows, with too large error bars. 

4 Results from local correlation tracking 

We also apply the local correlation tracking (hereafter LCT) 
technique, which was originally designed to remove seeing- 
induced distortion in sequences of solar images. This tech- 
nique examines the best match between the pixels of the 
consecutive images in a correlation window. The best match 
defines a vector of displacement, which, as a result gives 
a map of surface flow patterns. For further details of the 
method as well as for application see Svanda et al. (2006, 
2007) and our forthcoming paper for numerical tests on ar- 
tificial data by Svanda et al. (2008). 

The resulting surface flow map for the 1998 data is plot- 
ted in Fig.|2] indicating a complex flow pattern with conver- 
gent and also divergent flows, and maybe vortex-like fea- 
tures around the spots near 100°. The differential rotation 
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Fig. 2 The result of the LCT technique showing the surface flow pattern on UZ Librae from the time-series Doppler 
images in 1998. As a background the average surface temperature map is plotted. 




pattern can be derived from averaging the zonal flow com- 
ponents. The best-fit numerical result suggests a weak but 
uncertain anti-solar shear of —0.001 ± 0.003 which is at 
least not in a contradiction with the ACCORD results. The 
latitude shift - longitude diagram from LCT is overplotted 
on the corresponding ccf-map in Fig.[T]with triangles. Un- 
fortunately, the LCT technique cannot be applied for the two 
separated subsets in 2000. 

From averaging the meridional flow components we are 
able to investigate the meridional motions. The resulting 
plot is seen in Fig. [5] (bottom). We get small velocity values 
in the order of a few hundreds ms^^, and the comparison 
with the ACCORD result is meaningless. 

5 Summary 

- Applying our cross-correlation technique ACCORD for 
the time-series Dopper images of UZLib in 1998 and 
in 2000 we derived weak surface DR of anti-solar type, 
in agreement with former DR detections. However, no 
sign of coherent meridional flow was found. 

- LCT technique is applied first time for stellar data (see 
also the paper by Kovari et al. 2007a) and delineated a 
complex surface flow network for UZ Lib. 

- Averaging the zonal components of the LCT flow map 
along latitude resulted in an uncertain, weak anti-solar 
DR. Meridional flow components show no evidence for 
any coherent pattern, similarly to the ACCORD non- 
detection. 
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Fig. 3 Best correlated latitudinal shifts from ACCORD 
(top) and meridional flow components from LCT (bottom). 
No sign of coherent meridional flow can be seen. 
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